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ABSTRACT 

The ultra-long period Cepheids (ULPCs) are classical Cepheids with pulsation periods exceeding 
^ 80 days. The intrinsic brightness of ULPCs are ^ 1 to ^ 3 mag brighter than their shorter period 
counterparts. This makes them attractive in future distance scale work to derive distances beyond the 
limit set by the shorter period Cepheids. We have initiated a program to search for ULPCs in M31, 
using the single-band data taken from the Palomar Transient Factory, and identified eight possible 
candidates. In this work, we presented the U/-band follow-up observations of these eight candidates. 
Based on our U/-band light curves of these candidates and their locations in the color-magnitude 
diagram and the Period-Wesenheit diagram, we verify two candidates as being truly ULPCs. The six 
other candidates are most likely other kinds of long-period variables. With the two confirmed M31 
ULPCs, we tested the applicability of ULPCs in distance scale work by deriving the distance modulus 
of M31. It was found to be ij^m 31 ,ulpc = 24.30 ± 0.76 mag. The large error in the derived distance 
modulus, together with the large intrinsic dispersion of the Period-Wesenheit (PW) relation and the 
small number of ULPCs in a given host galaxy, means that the question of the suitability of ULPCs as 
standard candles is still open. Further work is needed to enlarge the sample of calibrating ULPCs and 
reduce the intrinsic dispersion of the PW relation before re-considering ULPCs as suitable distance 
indicators. 

Subject headings: Cepheids — variables stars — galaxies: individual (M31) — distance scale — stars: 
distances 


1. INTRODUCTION 

Cepheid variables span a class of pulsating stars within 
the instability strip in the Hertzsprung-Russell dia¬ 
gram. Their well-known period-luminosity (PL) rela- 
ti on (also known as the Le avitt law, as first presented 
in lLeavitt fc PickerindllQll) makes Cepheids good stan¬ 
dard candles. They constitute the first rung of the 
extra-galactic distance scale ladder. Previous stud¬ 
ies have shown that the Cepheid PL relation can be 
used to deter mine distances out to the order of 10- 
30 Mpc (e.g., iFreedman et al.l [2QQ1I: iRiess et ^ 1201 1[) 
using the Hubble Space Telescope (HST). When cou¬ 
pled with other secondary distance indicators, such as 
type la supernovae and the Tully-Fisher relation, one 
is able to determine distances to galaxies that are well 
within the Hubble flow. However, the cosmic distance 
ladder suffers fro m several uncert ainties at each lad¬ 
der rung. Thus, iBird et ahl (|2009[ ) have proposed the 
use of very luminous Cepheids, the so-called ultra-long- 
period Cepheids (hereafter ULPCs), to extend distance 
estimation beyond ~ 100 Mpc. Besides the potential 
to be used as distance indicators, ULPCs can also be 
used to constrain and enhance our understanding of stel¬ 
lar pulsation and evolution theories for inter mediate- to 
high-mass stars crossing the i nstability st r ip (jBird et ^ 
120091: iFiorentino et al.ll2012[) . iBird et ^ (|20Q9[ ) defined 
ULPCs to be fundamental-mode (FU) Cepheids with pe- 
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riods longer than ^ 80 days. 

Most studies of extragalactic distance determination 
use the Large Magellanic Cloud (LMC) as a distance 
anchor. Since the LMC has irregularity in its three- 
dimensional shape and low metallicity, many authors 
suggest instead using M31 as a st epping stone to cosmic- 
distance determinat ion (see, e.g.. iClementini et al.ll200l1: 
IVilardell et al.|[2QTQl and references therein). The merits 
of M31 include the following: it has a simple geometry; 
potential M31 distance-indicator stars are bright enough 
to be resolved; it is a local counterpart to the spiral galax¬ 
ies that a re used to determine t he extragalactic distance 
(see, e.g.. IFreedman et ~al]l2001[ ): and it is a local bench¬ 
mark to calibrate the Tully-Fisher relation. Despite its 
proximity, M31 has not been cons idered as a distance 
anchor (see, e.g., IRiess et al.l |2Q1 ill . because a distance 
anchor has to have a precise distance estimate and har¬ 
bor a fair amount of distance indicators such as Cepheids. 
M31 has not yet met these two criteria because (1) its 
current distance est imate has larger uncertainty (^ 4%, 
IVilardell et al.ll^QlQl) compared to othe r distance anchors 
(for e xample, LMC at ^ 2% level; see iPietrzvhski et al.l 
I2Q13D . and (2) there was no large sample of well stud¬ 
ied Cepheids in the literature. Nevertheless, r ecent 
studies of M31 eclipsing binaries (iLee et al.|[2Q14b) and 
Cepheids (|Firi fc Valls-Gabaudl 1^121: IRiess et al.ll2Q12l: 
iKodric et 2Q13L 120141) have demonstrated the poten¬ 
tial of establishing M31 as a distance anchor in the near 
future. Given the increasing importance of M31 in fu¬ 
ture distance scale work, it makes sense to search for 
and identify ULPCs in M31: such studies could provide 
a “one-step” calibration of the Hubble consta nt, similar 
to the role of NGC 4258 in IRiess et al.l (|2QllD . 

Ground-based CCD observations of Cepheids in M31 
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Figure 1. Footprints for some of the previous surveys, as summarized in Table [T] superimposed on a PTF image of M31. Note that one 
CCD (upper row, third from the left) is inoperable. 


originate either from dedicated surveys to search for 
Cepheids and (detached) eclipsing binaries, or as by¬ 
products from intense monitoring of M31 to detect 
micr o-lensing events. Examples of the former case in ¬ 
clud e iFreedmanfc Madore (1 1 9901) . iMagnier et al.l (|1997[ ) , 


and IVilardell et al. 


well as the DI¬ 


RECT project fsee lKalnznv et anil998L and subsequent 
papers in the series). For the micro-lensing exper i- 


2003) . 


120101) . 


ments, these include the AGAPE (lAnsari et al.l 
the Nainital Microlensing S urvey (iJoshi et ahl 
the POINT-A GAPE Survey (1^ et al. 1 12003) . the We- 
CA PP survey (|Fiiri et al.ll2Q06l) . and the POMME Sur¬ 
vey (|Fliri fc Valls-Gabandll2012[) . In addition, time-series 
observations of two M31 fields, using the Large Binocular 
Telescope, have also detected a n umber of s hort-period 
Cephe ids (|Clementini et al.ll2011[) . Recently. iRiess et al.l 
(|2012[ ) combined data from the POMME Survey and 
HST observations for 68 Cepheids, and derive a true dis- 
tance of 752 ± 27 kpc to the M31. iKodric et al.l (j2013D 
presented a catalog of M31 Cepheids, including funda¬ 
mental and first overtone Cepheids (as well as Type II 
Cepheids), based on the first year of the Panoramic Sur¬ 
vey Telescope and Rapid Response System (PSl, Pan- 
STARRSl) PAndromeda Survey. A summary of these 
surveys and projects is presented in Table [U and the 
coverage for some of the surveys is shown in Figure [TJ 
These studies are limited in that most of the observa¬ 
tions, with the exception of the POMME Survey and the 
PAndromeda Survey, cover only part of M31, and most 
of them concentrate on the disk (see Figure [1]). Further¬ 
more, no Cepheids with periods longer than ^ 80 days 


were detected in these studies. 

Therefore, we have initiated a program to search for 
ULPCs in M31 by using data from the Palomar T ransient 
Factory (PTF: ll^w et al.ll2QQ9l:iRan et al.ll2QQ9[) . This is 
because the PTF data can cover the entire disk of M31, as 
shown in Figure [H In addition, PTF observes M31 rou¬ 
tinely with cadence up to one day. These two conditions 
make PTF an ideal data set to search for t he ULPCs 
i n M3 1. Our search results w ere published in iLee et al.l 
(|2Q13L hereafter Paper I) and iLee et al.l (|2Q14aD . Using 
image subtraction techniques, eight ULPC candidates 
were identified. Figure [2] displays the postage stamp im¬ 
ages for these candidates based on PTF data. As men¬ 
tioned in Paper I, time series U/-band data is needed 
to further confirm or disprove the ULPC nature of these 
candidates. In this work, we report the U/-band follow¬ 
up observations for all these eight candidates. Section 
2 describes the follow-up observations, data reduction, 
and the calibration of the light curves. Analysis of these 
light curves and the results are presented in Section 3, 
followed by our conclusions in Section 4. 

2. THE FOLLOW-UP OBSERVATIONS 

The U/-band follow-up obs ervations were carr ied out 
with the P60 Telescope (P6Q: ICenko et allDOOGl) at the 
Palomar Observatory and the Lulin 1 m Telescope (LOT; 
iKinoshita et al.|[2QQ5[) at Lulin Observatory. These ob¬ 
servations began in 2012 October and ended in 2013 
November. P60 is equipped with a 2k x 2k CCD with 
a pixel scale of 0.379 arcsec per pixel, while LOT used 
a Ik X Ik CCD with a pixel scale of 0.512 arcsec per 
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Table 1 

Summary of Modern-day Cepheid Observations in M31. 


Project/Survey^ 

Telescope'^ 

Field Center^ 

Field-of-View 

Pix. Scale*^ 

Obs. Time Span® 

Ref.f 


Dedicated Surveys to Search for Cepheids (and Eclipsing Binaries) 

Freedman & Madore 

3.6m CFHT 

10 fields covering Baade’s Field I, III & IV 

2' X 3'g 


1985-1988 

1 

Magnier et al. 

2.5m INT 

9 fields along spiral arms 

12.5' X 12.5' 

0.367 

1993 (9) 

2 


1.3m MDM 

(same as above) 

10.9' X 10.9' 

0.637 

1993 (11) 

2 

Vilardell et al. 

2.5m INT 

q; = 00^44^^46^ S = +41^38'20" 

33.8' X 33.8' 

0.33 

1999-2003 (21) 

3, 4 

DIRECT: M31B 

Two Tel.^ 

a = 11°.20 S = +41°.59 

11' X 11' 

0.32 

1996-1997 (~43) 

5 

DIRECT: M31A 

Two Tel.^ 

a = 11°.34 5 = +41°.37 

11' X 11' 

0.32 

1996-1997 (~44) 

6 

DIRECT: M31C 

Two Tel.^ 

a = 11°.10 6 = +41°.42 

11' X 11' 

0.32 

1996-1997 ^53) 

7 

DIRECT: M31D 

Two Tel.h 

a = 11°.03 5 = +41°.27 

11' X 11' 

0.32 

1996-1997 Css) 

8 

DIRECT: M31F 

Two Tel.^ 

a = 10°.10 5 = +40°.72 

11' X 11' 

0.32 

1996-1997 (~51) 

9 

DIRECT: M31Y 

1.2m FLWO 

a = 10°.97 5 = +41°.69 

11' X 11' 

0.33 

1999-2000 (~25) 

10 

By-products from Micro-lensing Experiments 

AGAPE 

2.0m TBL 

6 + 1 fields centered at M31, 

14' X 10'‘ 

0.30 

1994-1996 (~ 69) 

11 



oriented along the main axis 





NMS 

1.0m ST 

Q = 00''43“38= 5 = +41°09.1' 

13' X 13' 

0.37 

1998-2002 (> 150) 

12, 13 

POINT-AGAPE 

2.5m INT 

a = 00’’43™ 10® S = +40°58'15.0" 

33' X 33' 

0.33 

1999-2001 C ISO) 

14 



a = 00'’44™00® 5 = +41°34'00.0" 

33' X 33' 

0.33 

1999-2001 C ISO) 

14 

WeCAPP 

1.2m CAO 

a = 00*'42™44®.3 S = +41°16'07.5" 

17.2' X 17.2' 

0.50 

2000-2001 

15 


0.8m WO 

mosaic CAO’s FOV with 4 pointings 

8.3' X 8.3' 

0.49 

1999-2008 

15 

POMME 

3.6m CFHT 

Q = 00*'43™50® S = +41°45'0" 

h-^ 

0 

X 

h-^ 

0 

0.187 

2004 (~ 50) 

16 



a = 00h41™50® 6 = +40°44'0" 

h-^ 

0 

X 

h-^ 

0 

0.187 

2005 C 50) 

16 

Other Time-series Observations 

Clementini et al. 

8.4m LBT 

a = 00*'48™13.11® 5 = +40°19'09.4" 

23' X 23' 

0.225 

2007 (~ 8) 

17 



a = 00''49™08.31® 5 = +42°16'09.4" 

23' X 23' 

0.225 

2007 C S) 

17 

PAndromeda 

1.8m PSl 

a = 00*'42™44.33® 5 = +41°16'07.5" 

~ 2.6°x ~ 2.6° 

0.258 

2010-2011 (~ 183) 

18 


^ AGAPE = Andromeda Gravitational Amplification Pixel Experiment; NMS = Nainital Microlensing Survey; POINT-AGAPE = Pixel-lensing Obser¬ 
vations with the Isaac Newton Telescope-Andromeda Galaxy Amplified Pixels Experiment; WeGAPP = Wendelstein Galar Alto Pixellensing Project; 
POMME = Pixel Observations of M31 with MEgacam 

^ GFHT = Ganada-France-Hawaii Telescope (Hawaii, USA); INT = Issac Newton Telescope (Spain); MDM = McGraw-Hill Telescope at the Michigan- 
Dartmouth-MIT (MDM) Observatory; FLWO = F. L. Whipple Observatory (Arizona, USA); TBL = Bernard Lyot Telescope (France); WO = Wendelstein 
Observatory (Germany); GAO = Galar Alto Observatory (Spain); ST = Sampurnanand Telescope (India); LBT = Large Binocular Telescope (Arizona, 
USA); PSl = Pan-STARRSl Telescope (Hawaii, USA) 

^ R.A. and decl. of field center is in J2000. 

^ GGD pixel scale in ’’/pixel. 

® Number in the parenthesis is the number of observing nights when available. 

^ Reference: 1. iFreedman &: Madord ([199^: 2. iMagnier et al.l (| 1997li: 3. iVilardell et al.l |2QQ^; 4. iVilardell et al.l (I2QQ7I) ; 5. iKaluznv et al.l (| 199^: 6. 
IStanek et al.l (119981) : 7. IStanek et al.l (1^9^ : 8. IKaluznv et al.l (119991) : 9. iMocheiska et all (^9^ : 10. iBonanos et al.l ^20031) : 11. lAnsari et al.l (l2Q04f) : 12. 
iJoshi et al.1 (l2QQ^; 13. iJoshi et al.l (I2QT^ : 14. lAn et al.l (l2QQii : 15. iFliri et al.l (l20Q^ : 16. iFliri fc Valls-GabaudI (l2QTli : 17. IGlementini et al.l (I2QT1) : 18. 
IKodric et al.l (120131) 

S The FOV for each of the 10 GGD fields. 

^ The two telescopes used in the DIREGT project are: 1.3m MDM and 1.2m FLWO. 

^ This is the total FOV; the FOV for individual fields is 4' x 4.5'. 



Figure 2. Postage-stamp images of the eight ULPC candidates reported in Paper 1. The upper row shows the PTF R-band images, while 
the lower row shows the PTF ^f-band images. Note that the R- and ^f-band images are from different epochs. All stamps have a FOV of 
30" X 30". The green circles with a radius of 5" indicate the locations of the candidates. 
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Figure 3. Seeing distributions in the V (upper panel) and / bands 
(lower panel) for the P60 and LOT images. For the P60 observa¬ 
tions, besides a pair of y/-band images, an additional V-band 
image was also taken on the nights between 2012 November and 
2013 March. 

pixel. Observations for both telescopes were obtained 
in queue mode when there was available observing time 
in suitable weather conditions. Imaging data from P60 
was reduced with a dedicated pipeline as described in 
iCenko et al.l (|2QQ6[) . For LOT images, subroutines in 
IRAF0 were used to reduce the data in the usual manner 
(bias and dark subtraction, and flat-fielded with mas¬ 
ter flat images). The LOT /-band images suffered from 
a fringing problem. To remove the fringe patterns, we 
median combined the /-band images for candidates with 
an internal ID 4-10470 and created a master fringe im¬ 
age. We selected the images for this candidate because 
it is located far away from the M31 bulge and disk (see 
Figure 1 in Paper I). Hence, the images do not have a 
strong background gradient. The master fringe image 
was smoothed by a “boxcar” averaging algorithm avail¬ 
able from IRAF, and then scaled and subtracted from all 
of the LOT /-band images. Astrometric refinement for 
both P60- and LOT-r ednced images we re performed us¬ 
ing astrometry.net (|Lang et al.ll2QlQl) . Figure [3] shows 
the seeing distributions of all the images taken from both 
telescopes. It can be seen that some images with large 
seeing could be affected by bad weather. After visual in¬ 
spection, images that were affected by bad weather were 
discarded. This left 568 and 386 images from P60 and 
LOT, respectively. 

Photometry based on point-spread function fitting 
(hereafter PSF photometry) for the candidates and 
suitable stars in each frame were obtained from 
IRAF/DAOPHOT subroutines. For each candidate, we first 
created a “master” catalog (with ^ 50 stars) that in¬ 
cluded: (1) the location of the ULPC candidate itself; 
(2) faint stars in the vicinity of the candidate star (se¬ 
lected based on the best seeing images); and (3) relatively 

^ IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

® See Paper I for the meaning of internal ID. 


Table 2 

U/-band Light Curves for the Candidates. 


Candidate 

Telescope 

Band 

MJD 

m 


8-0326 

LOT 

V 

56218.62510 

18.819 

0.031 

8-0326 

LOT 

V 

56224.71680 

18.931 

0.034 

8-0326 

LOT 

V 

56224.71403 

18.935 

0.033 

8-0326 

LOT 

V 

56235.64163 

19.026 

0.038 


Note. — (This table is available in its entirety in machine-readable 
and Virtual Observatory (VO) forms.) 


bright (Vlgs ^ 16.3mag. to ^ 19.9mag.) and isolated 
stars t aken from the M31 Local Group Survey catalog 
(LGS; iMassev et al.l |2QQ^ . These LGS stars were vi¬ 
sually inspected to ensure they were located within the 
P60 and LOT images, as well as away from any crowded 
regions in the images. About ^ 15 to ^ 25 LGS stars 
in the “master” catalog were used to construct the PSF 
model for each image by executing the IRAF subroutines 
DAOPHOT.PHOT, DAOPHOT.PSTSELECT and DAOPHOT.PSF. 
Sky coordinates for all the stars in the “master” cat¬ 
alogs were then converted to pixel coordinates and 
saved as a coordinate file. Finally, DAOPHOT.PHOT and 
DAOPHOT. ALLSTAR were run with the input coordinate 
file and PSF model to obtain the instrumental PSF pho¬ 
tometry for the stars in the “master” catalog. 

2.1. Light Curves from Differential Photometry 

The R/-band light curves for the eight candidates 
were construct ed using differenti al photometry tech¬ 
niques given in iBroeg et al.l (|2QQ5f ) . Instead of using a 
single star as a comparison star (cs) to obt ain differen¬ 
tial ph otometry, the technique presented in iBroeg et al.l 
(|2QQ5[) used a number of stars (either all stars in the im¬ 
ages or a subset of suitable constant stars) to construct 
an artificial comparison star: 

Am = rUc—< m >cs = mg — Wjmf ^ (1) 

where me was the instrumental magnitudes for our candi¬ 
dates in either the V oi I band, and Wi were the weights 
for individual comparison stars such that ^ • Wi = 1. At 
the first stage, we constructed differential photometry 
light curves for a pair of LGS stars in our “master” cata¬ 
logs, as mentioned in the previous section, by taking one 
of them as a comparison star. We then removed those 
LGS stars that exhibited a large scatter in the differ¬ 
ential light curves. This left a subset of good LGS stars 
before equation (1) was applied to each of our candidates. 
Since the calibrated R/-band magnitudes were available 
for these comparison stars from the LGS catalog, light 
curves constructed from equation (1) were calibrated via 
the following equation: 

^calibrate = Ato + ^ , (2) 

i 

where Wi were the same weights as in equation (1), and 
m^S’LGS R/-band magnitudes from the LGS cat¬ 

alog. The calibrated R/-band light curves for the eight 
candidates are given in Table O and displayed in panel 
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(a) and (b) of Figures |1][71 

3. ANALYSIS AND RESULTS 
3.1. Period Search and Folded Light Curves 

The new U/-band light curves suggested that periods 
for the candidates need to be revised from Paper L For 
each pair of the U/-band light curves, we first subtracted 
the arithmetic mean of the light curve, followed by nor¬ 
malization of the light curve using the peak-to-peak am¬ 
plitude. We then applied the well-documented Lomb- 
Scargle algorithm to search for periodic signals by com¬ 
bining these mean-subtracted and normalized pair ofVI- 
band light curves (as shown in panel (c) of Figures |4]17j) 
to each of the candidates. We only searched for periodic 
signals around ^ ±30 or ^ ±50 days from the periods 
reported in Paper I (with the exception of candidate 8- 
0272, for which a shorter period range was searched for 
periodic signals). The periods corresponding to the high¬ 
est peaks in the Lomb-Scargle periodogram were taken 
as the periods for these candidates. They are indicated 
by a vertical dashed line in panel (d) of Figures HE as 
well as listed in the second column of Table [H We also 
ran a Monte-Carlo simulation to quantify the errors on 
the corresponding periodic signals. Assuming the photo¬ 
metric errors are Gaussian, about 10,000 simulated light 
curves were generated using the Box-Muller method for 
each of the candidates. The Lomb-Scargle algorithm was 
then applied to these light curves in a similar fashion to 
the combined U/-band light curves. The errors on the 
derived periods were then estimated as the rms of the 
corresponding peaks in the power spectra. 

Using the periods given in Tabled we folded the VI- 
band light curves for these eight candidates. The folded 
U/-band light curves were then fitted wit h a low-order 
Fouri er decomposition (for example, see iSimon fc Led 
|1981[) in the following form: 

i=2,3,4 

m(0) = mo ± ^ Ai cos(2i7r0 ± (3) 

where 0 G [0,1] were phases of the pulsating cycles. 
The fitted light curves, shown as black curves in panels 
(e) and (f) of Figures |4]171 were used to derive the VI- 
band mean magnitudes (by converting the fitted light 
curves to intensity, taking the average, and converting 
back to magnitude), as well as the extinction-free We- 
senheit func tion W. In order to be consistent with pre¬ 
vious works (|Bird et al.ll2QQ9l: iFiorenUno et al.ll2Q12[) . we 
adopted the same Wese nheit function in the form of 
W = I — l.bb(V — I) from FUdalski et al.l (|1999al) . The co¬ 
efficient of 1.55 in th e Wesenheit fu n ction is based on the 
extinction law from iSchlegel et al.l (|1998[) . The derived 
mean magnitudes and the values of Wesenheit function 
are listed in Table [3l 

The shape of the ULPC U/-band light curves are 
similar to the c lassical Cepheids (|Bird et al.l I2QQ9I: 
iNgeow et al.l I2Q13F for which light curves exhibit the 
characteristic “saw-tooth” shape with a quick rise to 
maximum light followed by a slow decline to minimum 
light. Based on our U/-band light curves, only two of 
the candidates (8-0326 and 8-1498) exhibit Cepheid-like 
light curves, as shown in Figures S] and [5l respectively. 


Table 3 

Revised Periods and Mean Magnitudes for the Candidates. 


Candidate 

Period (days) 

< U > 

< / > 

W 


M31 ULPC Candidates 

8-0326 

74.427 ±0.120 

18.684 

17.256 

15.043 

8-1498 

83.181 ±0.178 

18.856 

17.783 

16.120 

M31 non-ULPC Candidates 

8-1176 

125.313 ±4.832 

19.095 

17.367 

14.689 

7-1326 

132.608 ± 4.543 

19.357 

17.356 

14.254 

8-0272 

169.262 ±4.461 

20.464 

18.176 

14.630 

4-1047 

216.169 ±0.907 

18.625 

17.548 

15.879 

9-0530 

221.779 ± 17.695 

20.218 

19.092 

17.347 

8-1180 

222.124 ±5.959 

19.780 

17.246 

13.318 


The U/-band light curves for the other six candidates 
are given in Figures [6] and [71 These curves do not dis¬ 
play light curves that are expected for Cepheids. In the 
following sub-sections, the two candidates with Cepheid- 
like light curves are referred to as “M31 ULPC candi¬ 
dates,” while the other six candidates will be referred to 
as “M31 non-ULPC candidates.” Further confirmation 
or falsification of their ULPC nature will be detailed in 
the next sub-section. We note that the two “M31 ULPC 
candidates” have well-determined periods (with errors of 
the order of ^ 0.2%), while the periods for the “M31 non- 
ULPC candidates,” with the exception of the candidate 
with an internal ID of 4-1047, have a much larger error 
(from ^ 2.6% to ^ 8%). Furthermore, as mentioned in 
Paper I, these two “M31 ULPC candidates” have pre¬ 
viously bee n classified as Cephe ids in the literature (for 
example, in iMagnier et al.iri997l but they have incorrect 
periods). 

The two M31 ULPC candidates have periods of ^ 74.4 
and ^ 83.2 days, r espec tively. According to the defini¬ 
tion by [Bird et al.l (|2Q09[ ). a ULPC should have a period 
longer than 80 days. Hence, the ^74.4 days ULPC found 
in this work may not be classified as ULPC. However, an 
80 day period cut is rather arbitrary, and we prefer to 
identify the ULPCs based on the discussion given in Sec¬ 
tion 3.2. Furthermore, initial period detection searches 
for this candidate gave results in excess of 80 days in 
Paper I, and was refined to a shorter period with other 
period-search algorithms and U/-band light curve data. 
Therefore, we retained this candidate in the sample of 
M31 ULPCs. 

3.2. Verification of the ULPC Candidates 

Besides the shapes of the U/-band light curves, we 
verify/falsify the ULPC nature of our eight candidates 
by comparing them to the known ULPCs (given in 
iFiorentino et al.l |2012[) in the color-magnitude diagram 
(CMD) and the Period-Wesenheit diagram (PWD), col¬ 
lectively presented in Figure [8l Using the CMD and 
PWD to identify Cepheids has b een done, for exam - 
ple, in lUdalski et al.l (|1999bD and iKodric et al.l (|2013D . 
Only those Cepheids identified with pe riods longer than 
^ 80 days will be re-cl assified as ULPC (|Bird et al.ll2009l: 
IFiorentino et al .112 01 21) . When making such comparisons, 
we have adopted the distance modulus of M31 as rm.s i = 
24.46 mag, recommended bv ide Griis fc Bonol (|2014f ). to 
convert the observed magnitudes of these candidates to 
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Figure 4. (a) Observed V-band light curve for candidate 8-0326, where the green open circles and red filled squares are for the P60 and 
LOT data, respectively, (b) Observed /-band light curve for the same candidate, (c) the combined VI band light curve that was used 
in period search; see the text (Section 3) for details, (d) The Lomb-Scargle periodogram, at which the corresponding period at the peak 
of the periodogram (indicates by a vertical dashed line), is the adopted period listed in Table O (e) The V-band folded light curve using 
the period given in Table [3] (f) The /-band folded light curve for the same candidate. The black curves in (e) and (f) represent the fitted 
light curves based on low-order Fourier decomposition. 
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Figure 5. Same as Figure lU but for candidate 8-1498. 

the absolute magnitudes. We also included the LMC 
FU Cepheids and long-period variable s (LPV) based 
on the OG LE-III database , taken from iSoszvnski et al.l 
(|2QQ8[ ) and ISoszvnski et al.l (j2QQ9f ) respectively, in these 
diagrams. Again, we adopted jiLMC = 18.49 mag 
(|de Grijs et al .ll2()T4 ) when converting the observed mag¬ 
nitudes to absolute magnitudes. We enlarged Figure [8] by 
only showing the known ULPGs and our eight candidates 
in Figure [9l Figures [HI and [9] reveal that the two “M31 
ULPG candidates” are indeed ULPG, as both of them fall 
in the parameter spaces defined by the known ULPGs. 
The remaining six “M31 non-ULPG candidates” are fal¬ 
sified as ULPG based on their location in the GMD and 
PWD diagrams, even though a few of them may be con¬ 
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sistent with ULPGs in either the GMD or PWD diagram, 
but not in both. For example, candidates 4-1047 and 9- 
0530 fall within the ULPG instability strip, but they are 
outside the ±3cr boundary of the Period-Wesenheit (PW) 
relation. Similarly, candidates 7-1326, 8-0272 and 8-1176 
are located within the ±3cr boundary of the PW relation, 
however they are too red in the GMD to be considered 
as ULPG. Since the main goal of this work is to verify 
the ULPG nature of our candidates in M31, determining 
the nature of these six non-ULPG candidates is beyond 
the scope of this paper, and they will not be discussed 
further. 

It is expected that a galaxy would have small num¬ 
ber of ULPGs for two reasons. Firstly, ULPGs are 
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Figure 6. Same as Figure lU but for candidate 4-1047, 8-0272 & 9-0530. 
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Figure 7. Same as Figure 01 but for candidate 7-1326, 8-1176 & 8-1180. We found that an additional zero-point offset of +0.078, +0.089 
and +0.078 is needed to be added to P60 data for these candidates with internal ID 7-1326, 8-1176 and 8-1180, respectively. 
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Figure 8. Comparison of our eight candidates to the known ULPCs listed in IFiorentino et "HI (120121) . as well as the LMC Cepheids and 
long-period variables (LPV) from the OGLE-III database, in the observed CMD (left panel) and PWD (right panel). We included LPVs 
to aid in falsifying the ULPC nature of our c andidates. The LPVs were classified as small amplitude red giants (OSARC), semi-regular 
variables (SRV), and Mira by the OGLE team (|Soszvnski et al1l2009l ). Note that the saturation limit in OGLE-III data is around 13 mag 
to ~ 13.5 mag. This translates to My ~ —5.0 mag to ~ —5.5 mag. 
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Figure 9. Observed GMD (left panel) and PWD (right panel) only showing the ULPGs. Since the predicted bound ary of the ULPG 
insta bility strip from theoretical calculations is not yet available, we used the color range available from known ULPCs (|Fiorentino et all 
120121) . 0.6 < (V-I) < 1.6, to represent the (crude) boundary of the i nstability strip (shown as vertical lines in the left panel). The 
solid line in the right panel is the Period-Wesenheit relation adopted from [FTorentino et al and the dashed lines represent the ±3(7 

boundary (where a = 0.34 mag). 


Table 4 

M31 Cepheids with Periods Greater than ~ 75 days but Less than 80 days. 


Name or ID 

Period (days) 

a(J2000) 

5(J2000) 

Reference 

vn.4.2.678 

PSO J010.5806±40.8319 

78.00 

74.79 

11.06420 

10.58066 

41.56927 

40.83192 

Riess et al. (2012) 

Kodric et al. (2013. onlv for Geoheids classified as “FU”) 


intermediate- to high- mass stars with masses between 
^ 15 Mc^ ang ^ 20M© (|Bird et al.ll2QQ9l:lFiQrentinQ et al.l 
I2Q12F Hence, their number should be lower than their 
shorter period counterparts that have lower mass. Sec¬ 


ond, these intermediate- to high-mass stars need to cross 
the instability strip in order to pulsate and the time that 
they spend inside the instability strip is relatively short 
as compared to their evolutionary age. Based on the 
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list compiled in iFiorentino et al.l (|2Q12[) , the number of 
ULPCs in a given galaxy ranges from one (NGC 6822), 
two (M81 and I Zw 18), three (SMC ang NGC 300), four 
(LMC, NGC 1309 and NGC 3021), five (NGC 55), and 
nine (NGC 3370). Then, it is not a surprise that only 
two ULPCs were found in M31. Another reason for de¬ 
tecting a small number of ULPCs in M31 could be due 
to characteristics of PTF data. For example, the one 
CCD that is out of commission covers a portion of M31’s 
disk (see Figure [T]), and the stars in M31’s bulge region 
are not resolvable in PTF images. Further discussion on 
the effect of PTF data can be found in Paper I, and will 
not be repeated here. A few additional ULPCs might be 
discovered from the three years of the Pan-STARRSl’s 
PAndromeda Survey (M. Kodric et al. 2014, private com¬ 
munication). On the other hand, if the period cut of 
^ 80 days that defines ULPC is lowered to ^ 75 days, 
then there are a few more very long-period Cepheids 
in M31 that could be classified as ULPCs. These are 
listed in Table IH They were not reported in Paper I be¬ 
cause their periods are lower than the search criterion of 
P > 80 days. Furthermore, there is no UJ-band photo¬ 
metric data available for these Cepheids currently. 

3.3. Distance Scale Application 

In the previous sub-section, we adopted a distance 
modulus for M31 in order to verify the ULPC nature of 
our candidates. In this sub-section, we reverse the prob¬ 
lem by using the two confirmed ULPCs to determine the 
distance modulus of M31, i.e., fi = W — Mw The abso¬ 
lute magnitude Mw of the two ULPCs c an be determined 
from t he latest PW relation as given in IFiorentino et ^ 
(j2Q12f ): Mw = —2.661og(P) — 3.68, with a dispersion of 
(7 = 0.34. Since the PW relation is a statistical relation, 
the dispersion of this relation (a) will dominate the er¬ 
ror term of the derived distance modulus for individual 
ULPCs. Hence, we have adopted the dispersion of the 
PW relation as the error in the calculated distance mod¬ 
uli for the two M31 ULPCs: /i(8-0326) = 23.70T0.34 mag 
and /i(8-1498) = 24.91 ± 0.34 mag. Since the weights for 
these two distance moduli are the same, taking an aver¬ 
age of them reduces to the case of an unweighted mean. 
This procedure yielded i^msi^ulpc = 24.30 ± 0.76 mag, 
where the error on the averaged distance modulus is 
calculated usin g the small number statistics given in 
iKeepind (jl962L p. 202). The large error on the deter¬ 
mined distance modulus is mainly due to the combina¬ 
tion of small number statistics (as there are only two 
ULPCs found in M31) and the large dispersion of the 
PW relation for ULPCs. 

G. Fiorentino (2014, private communication) suggested 
using the PW relation derived from ULPCs in metal-rich 
galaxies (hereafter metal-rich ULPCs) . However, the 
sampl e of metal-rich ULPCs given in IFiorentino et al.l 
(j2012f ). i.e. those with 12\og[0/H] > 8.4 dex, only oc¬ 
cupied a narrow range in log(P). This causes the derived 
PW relation to be unreliable. Instead, we have averaged 
the absolute Wesenheit magnitudes for this sample of 
metal-rich ULPCs, yielding Mw = —8.84 mag. Using 
this absolute Wesenheit magnitude, the derived distance 
modulus to M31 with our two ULPCs is iaM 3 i,ULPC = 
24.42 ± 0.68 mag with a (unbiased) standard deviation 
of 0.76. The value of this distance m odulus is almost 
identical to the value recommended in Ide Criis fc Bonol 



_J_^^^^^^^^^^^^^^^^^^^_L_ 

1.9 2 2.1 2.2 2.3 

log(P) 


Figure 10. Top: comparison of th e tw o Period-Wesenheit rela - 
tions adopted from IBird et al ] (|2009li and IFiorentino et al.l (|2012l f ; 
Bottom: difference in the derived distance modulus from these 
two Period-Wesenheit relations as a function of the ULPC periods. 

(IMl . Note that the san iple of metal-poor ULPC s (12 + 
\og[0/H] < 8.4 dex) in IFiorentino et a 1 (1201211 is the 
same as the sample given in iBird et al.l ( 2QQ9D . The de¬ 
rived distance modulus is I^msi^ulpc = 24.74 + 0 . 68 ma g 
when using the PW relation given in iBird et al.l (|2QQ9[ ). 

The di stance moduli derived using the PW re l ations 
given in IFiorentino et al.l (|2Q12D and iBird et all (|2QQ9f ) 
differed by 0.44 mag, albeit with large uncertainties. 
This reflects one of the current problems in using ULPCs 
to derive distances: the calibration of the ULPC PW re¬ 
lations is still uncertain with many discrep ant results. 
The s lope of the PW relation derived in IBird et ^ 
(I 200 I) . based on 18 ULPCs in metal-poor host galaxies, 
is essentially zero but with a larg e error (—0.05+0. 5 4). In 
contrast, the PW slope given in IFiorentino et aP (|2012f ) 
is —2.66 (though no error on the slope is quoted in their 
paper). The intercepts of these two Period-Wesenheit 
relations even displayed a larger difference: —9.06 + 1.12 
versus —3.68. The top panel in Figure [10] compares these 
two PW relations, while the bottom panel of the same 
figure shows the difference of the derived distance modu¬ 
lus based on these two PW relations: this can be as large 
as ^ 0.5 mag. This disagreement between the two PW 
relations suggests more work is needed in the future to 
properly calibrate the ULPC PW relati on. Nevertheless, 
the di stance modulus derived from the IFiorentino et al l 
(j2012f ) PW relation will be adopted in this work, be¬ 
cause this Period-Wesenheit relation is derived from a 
larger sample of ULPCs than can be found in the lit¬ 
erature. Figure [TT] shows a comparison of the adopted 
distance modulus (based on the two M31 ULPCs) to 
other distance moduli given in literature. Our calcu¬ 
lated distance modulus is consistent with some of re¬ 
cent determinations, including the recomm ended value 
of 24.46 ± 0.10 mag (|de Criis fc Bonol[2014j) . 

We note that the Wesenheit magnitudes and the dis¬ 
tance moduli for these two ULPCs differ by ^ 1 mag. 
This is not due to the “depth effect,” as one ULPC is 
located at the near-side and another ULPC is located at 
the far-side of M31 (see Figure [HI. Assuming the size 
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Figure 11. Comparison of the M31 distance modulus based on 
the two ULPCs in this work {iimzi,ulpc = 24.30 ± 0.76 mag, 
filled circle) to other distance moduli determined using Cepheids 
(as compiled in ide Griis &: Bonoll20ld . open triangles). Note that 
some of the distance moduli, especially those in early years, did not 
include an error estimation. The horizontal dash ed line represents 
the re commended distance modulus of 24.46 mag (Ide Griis &: Bonol 

EoID. 



18 12082 15965 19847 23730 27J 


Figure 12. Locations of the two ULPCs in M31. The M31 image, 
with a size of 60' x 60', was downloaded from the Digitized Sky 
Surveys (DSS) archive. 

of M31 is about 40 kpc, the maximum difference in dis¬ 
tance modulus for two stars located at the two extreme 
edges is ^ 0.11 mag. This is ^ 10x smaller than the 
observed difference. Instead, the observed ^ 1 mag dif¬ 
ference in Wesenheit magnitudes or distance moduli for 
these two ULPCs is consistent with the spread of Mw at 
log(P) ^ 1.9, as shown in right panel of Figure [8] (or Fig¬ 
ure 2 in iFiorentino et ani2012[) . and the expected spread 
of the PW relation (4 x a = 1.36 mag). We believe this 
is due to their relative locations within the instability 
strip on the CMD, depending on their evolutionary sta¬ 


tus H Figure [T3l compares the lo cations of the ULPCs i n 
individual galaxies, as listed in iFiorentino et aH (|2012[) . 
to the locations of the two M31 ULPCs in CMD and 
PWD. This shows that ULPCs in some galaxies (such as 
NGC 3370 and NGC 3021) exhibit a similar spread in 
Wesenheit magnitudes as in the case of the M31 ULPGs. 
The spread of Wesenheit magnitudes, at a given period, 
can be expressed as AMw = AM/ — 1.55A(U — /), 
where AM/ and A(U — I) represent the range of mag¬ 
nitudes and colors, respectively, for ULPCs in the same 
galaxy. Therefore, assuming AM/ ^ 0 and a difference 
ofA(U — /)^0.5 mag could translate to a difference of 
AMw ^ 0.8 which is close to the observed spread 
of M31 ULPCsff This implies that at a given period, the 
Wesenheit magnitudes or distance moduli for the ULPCs 
can differ by as much as by ^ 1 magnitude. This calls 
into question the use of ULPCs as standard candles. 

4. CONCLUSION 

In this work, we have presented the U/-band follow-up 
observations of eight ULPC candidates in M31. Based on 
their light curve shapes and their locations on the CMD 
and PWD, we verified that two of the candidates are 
indeed ULPCs, and the remaining six candidates are not 
ULPCs (and possibly belong to the LPV class). These 
six non-ULPC candidates also showed that other types of 
LPV could be mis-classified as ULPCs (and vice versa) if 
they have not had ap propriate follow-up observations in 
the U/-bands (see also iNgeow et al.ir2Q13[) . We then used 
these two confirmed ULPCs to test their applicability in 
distance scale work by deriving the distance modulus to 
M31. Our derived distance modulus is consistent with 
the recommended value given in the literature, but with 
a large error of 0.76 mag. We have demonstrated three 
problems when using ULPCs as a distance indicator. 

1. Small number statistics: the number of ULPCs is 
expected to be small in a host galaxy; hence, the 
ULPCs do not sample the instability strip well. 

2. The discrepancy between PW relations in the lit- 
eratnre: the tw o av ailable PW rela t ions from 
iBird et all (|2QQ9[ ) and IFiorentino et al.l (|2Q12[) are 
in disagreement. As a result, the derived distance 
moduli can differ by as much as ^ 0.5 mag. 

3. The large dispersion of the PW relation: the cur¬ 
rent calibration of PW relation, defined as VF = 
I — 1.55(U — /), for ULPCs still possesses a large 
dispersion compared to their shorter period coun- 
terparts (which i s ^ 5 x smaller as shown in 
iNgeow et al.ll2009l) . Furthermore, at a given pe¬ 
riod, the back-to-back scatter in the PWD can be 
as large as ^ 1 mag. 

Combining these three problems with currently available 
PW relations results in a less precise and accurate dis¬ 
tance modulus. As mentioned in Fiorentino’s oral pre- 

^ The evolutionary status of ULPCs could be affected by various 
physical parameters, such as age, metallicity, mass-loss, and etc. 
However, a detailed investigation of their evolutionary status is 
beyond the scope of this paper. 

® To be more specific, using the values give in Table [S] we found 
that AI = —0.527 mag and A(U — I) = 0.355 mag for the two 
M31 ULPCs yielding AMw ~ 1-08 mag as observed. 







































12 


Ngeow et al. 



0.6 0.8 1 1.2 1.4 0.6 0.8 1 1.2 1.4 0.6 0.8 1 1.2 1.4 


(v-i) 



Figure 13. Compari son of the locations of the two M31 ULPCs (red filled circles) with ULPCs in other galaxies (open squares, taken from 
IFiorentino et al1l2012l ') in the CMD (left panel) and PWD (right panel). We excluded the ULPC in NGC 6822, as this gal axy has only one 
ULPC . Note that sub-figures in each panel are ordered according to the metallicity of the host galaxy, which is given in IFiorentino et ahl 

(|20I2). 


sentatior0, ULPCs are not yet ready to be used as a 
standard candle. Larger samples of ULPCs with well 
calibrated distance host galaxy distances are needed in 
order to derive a better PW relation (with smaller a) 
before using ULPCs in “one-step” determinations of the 
Hubble constant in the future. 

Regardless of their potential or problems in future 
distance scale applications, ULPCs represent a unique 
probe for studies of stellar pulsation and evolution be¬ 
cause they occupy the upper part of the instability strip 
that has attracted little attention to date. Evolution¬ 
ary tracks based on stellar evolution models appropri - 
ate for ULPCs have been ex plored in iBird et al.l (j2QQ9D 
and IFiorentino et al.l (j2Q12[ ) and will not be repeated 
here. In Figure [CT we showed the isochrones adopted 
from iBertelli et al.l ~ (j2QQ9[lF1 on the CMD, suggesting that 
ULPCs are young objects as expected. Nevertheless, de¬ 
tailed investigation of the evolutionary status for ULPCs 
with theoretical evolutionary tracks and isochrones is be¬ 
yond the scope of this paper and will be addressed in 
future work. 
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